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Cohesive laws are stress-strain curves used in finite element calculations to describe the debonding 
of interfaces such as grain boundaries. It would be convenient to describe grain boundary cohesive 
laws as a function of the parameters needed to describe the grain boundary geometry; two parameters 
in 2D and 5 parameters in 3D. However, we find that the cohesive law is not a smooth function of 
these parameters. In fact, it is discontinuous at geometries for which the two grains have repeat 
distances that are rational with respect to one another. Using atomistic simulations, we extract 
grain boundary energies and cohesive laws of grain boundary fracture in 2D with a Lennard- Jones 
potential for all possible geometries which can be simulated within periodic boundary conditions 
with a maximum box size. We introduce a model where grain boundaries are represented as high 
symmetry boundaries decorated by extra dislocations. Using it, we develop a functional form for 
the symmetric grain boundary energies, which have cusps at all high symmetry angles. We also 
find the asymptotic form of the fracture toughness near the discontinuities at high symmetry grain 
boundaries using our dislocation decoration model. 

PACS numbers: 62.20.Mk, 61.72.Mm, 31.15.Qg 



I. INTRODUCTION 

In materials such as silicon and various aluminum al- 
loys, cracks initiate and propagate along the interfaces 
between polycrystals known as grain boundaries. When 
the cracks initiate at the site of these microscopic de- 
fects, the macroscopic fracture strength of the material 
is dependent on the microscopic structure of the grain 
boundaries. The debonding of an interface such as a 
grain boundary is described by a cohesive law, giving the 
displacement across the interface as a function of stress 
(figure [1]). 

Cohesive laws are used by finite element cohesive zone 
models, which simulate fracture initiation at interfaces^. 
It has been shown that the shape and scale of the co- 
hesive law has a large effect on the outcome of the fi- 
nite element simulationii^. However, the CZM studies of 
grain boundary fracture have used cohesive laws that are 
guessed, chosen for numerical convergence, and do not 
take into consideration the effect of varying grain bound- 
ary geometries within the material - the same cohesive 
law is often used throughout the material despite the fact 
that in a real material, grain boundaries of varying ge- 
ometries must occur— i^i^. 

It would be useful to find a formula for the cohesive 
laws of the grain boundaries of a given material as a func- 
tion of geometry, for input into finite element simulations. 
The geometry of a 3D grain boundary depends on 5 pa- 
rameters that describe the orientations of the two grains. 
In addition, there are three different modes of fracture 
(normal to the crack plane, shear in a direction parallel 
to the crack line, or shear in a direction perpendicular 
to the crack line) to explore^^ as well as dependencies on 
temperature, impurities at the interface, and emission of 
dislocations to consider. 

Thus far, no systematic study of grain boundary cohe- 



sive laws as a function of geometry has been done with 
molecular dynamics or experiment. There is difficulty in 
measuring cohesive laws experimentally because it is dif- 
ficult to isolate and measure the displacements on either 
side of the grain boundary. 
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FIG. 1: An Example of a Cohesive Law. The stress 
vs. strain curve that describes the debonding of a 2D grain 
boundary with tilt angles 33.418° and 26.58°. Our measure- 
ments include the elastic response of the perfect crystal on 
either side of the grain boundary. In the results shown above, 
we have subtracted off the elastic response of the bulk in order 
to isolate the elastic response of the interface. This is de- 
scribed in detail in^ where cohesive laws of grain boundaries 
are used in finite element simulations of polycrystals. 

Previous atomistic studies of the mechanical response 
of grain boundaries have concentrated on a small number 
of symmetric grain boundaries in 3d5,7,8,9 ,io,ii ^ Explor- 
ing the complete picture of 3D grain boundary cohesive 
laws involves exploring a 5-dimensional space with three 
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modes of fracture. Because of this difficulty, we have 
taken a step back. We seek to systematically explore the 
cohesive laws for mode I fracture for all possible grain 
boundary geometries in 2D that can be simulated in peri- 
odic boundary conditions for a particular size and strain. 
We will initially be focusing on symmetric grain bound- 
aries and then expanding the picture to look at asym- 
metric 2D grain boundaries. One use of such data would 
be for finite element simulations of polycrystals^, as de- 
scribed above. For this purpose, it would be nice to find 
a functional form which describes the fracture strength 
of the grain boundary as a function of geometry. 

We find, however, that the fracture strength as a func- 
tion of tilt angles is discontinuous everywhere, with par- 
ticularly large jumps at special, high symmetry grain 
boundaries composed of a simple arrangement of struc- 
tural units with a low repeat distance. These special 
boundaries are also associated with cusps in the grain 
boundary energy. We will describe dependence of the en- 
ergy and the fracture strength near these special bound- 
aries by treating them as perfect, albeit complex, crystals 
with added dislocations to break the symmetry. 

Most studies use the coincidence site modeU^ to 
construct and classify special grain boundaries^iii^iiSiii^ 
while a few mention the significance of structural 
unitsiiSiii^iiiii^. The only 2D study we are aware of, the 
"bubble-raft" model, observes the structural units of sev- 
eral of the special, high angle grain boundaries for the tri- 
angular lattice and how patterns of structural units com- 
bine to create vicinal geometriesi^. Other 3D studies dis- 
cuss how grain boundaries with geometries close to those 
of special grain boundaries can be decomposed into the 
structural unit of the special grain boundary with added 
flawsi^ii^. We develop a systematic way of finding high 
symmetry geometries and show that the combinations of 
patterns of structural units at vicinal grain boundaries 
are key to understanding the dependence of energy and 
fracture strength on geometry. 

Sansoz and Molinari find the grain boundary energy by 
allowing the grains to relax together from an initial sepa- 
ration of a few angstroms^ while others perform a conju- 
gate gradients search^. We use a systematic method for 
explicitly imposing a relative shift between the grain and 
using atomistic relaxation for finding the global energy 
minimum similar to that used in^ii^. It is well estab- 
lished that there exist cusps in the grain boundary energy 
for special high angle grain boundarie9^i^ i^°i^^i^^'^^i^^'^^ . 
Recent studies of grain boundary constitutive properties 
focus on the response to sheari^ or compare shear and 
tension^iii. Warner et al. claim that the tensile response 
does not depend on the geometry of the grain boundary^. 
Others have seen jumps in the tensile fracture strength 
at high angle grain boundaries gcomctries^ii^. 

We will explain the cusps in energy and jumps in frac- 
ture strength by drawing an analogy between perfect 
crystals and high symmetry, high angle grain boundaries. 
The dislocation model of low angle grain boundaries gives 
a 9 log form for the grain boundary energy. Because vic- 



inal grain boundaries can be thought of as high symmetry 
grain boundaries with added flaws, the energy cusps at 
special grain boundaries have the same 0\og6 form. A 
similar argument will apply to the peak fracture stress. 
Just as adding dislocations to the perfect crystal adds 
a nucleation point for fracture, and therefore a discon- 
tinuity (jump down) in fracture strength, adding a flaw 
to a high symmetry, high angle grain boundary abruptly 
changes the local fracture strength by adding a potential 
nucleation site for fracture. 



II. PROCEDURE FOR CALCULATING GRAIN 
BOUNDARY ENERGY AND COHESIVE 
STRENGTH 

We calculate energies and cohesive strengths for grain 
boundaries using the DigitalMaterial^ package to per- 
form atomistic simulations. The potential we are using 
for 2D simulations is the Lennard Jones potential with a 
smooth, fourth order cutoff between 2.41 and 2.7 A. The 
ground state is a triangular lattice with a lattice constant 
of 1.11 A. 




FIG. 2: An Example of a Grain Boundary. The tilt 
angles are given by 6i and 02 and the repeat distances are 
Di = |(ai,fei)| and D2 — \{a2,b2)\. We have assumed the 
convention where 0° indicates an orientation with Miller in- 
dices (2, 1) for the vertical face of the grain boundary. 

We set up the grain boundary by initializing two rect- 
angular grains with the given rotations that define the 
grain boundary geometry we wish to measure. An exam- 
ple of the configuration of atoms that make up a grain 
boundary simulation is shown in figure [21 We use peri- 
odic boundary conditions in the y-direction in order to 
avoid edge effects which could potentially cause cracks to 
nucleate at the intersection with the surface^i. This has 
the disadvantage of only allowing geometries that have fi- 
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nite repeat distances. A method for finding these geome- 
tries is discussed in detail in section IIIIBI If necessary, 
both grains are strained shghtly by equal and opposite 
amounts in order to have both grains fit into a periodic 
box. We use a constrained layer of atoms to impose fixed 
boundary conditions in the x-direction. These are repre- 
sented by the darker atoms in figure [21 The constrained 
layer of atoms has a width equal to twice the cutoff dis- 
tance of the potential to eliminate surface effects from 
the free atoms. 

Besides the tilt angles, there are other factors to con- 
sider in constructing the grain boundary geometries. For 
commensurate grain boundaries, there will be an ideal 
relative displacement in the direction parallel to the grain 
boundary (the y-direction in figure [2]). This ideal dis- 
placement will correspond to the lowest energy and thus 
is the most natural configuration for a given pair of tilt 
angles 1^ 

We find the displacement along the boundary which 
gives the lowest energy by initializing the two grains 
with a small displacement in the x-direction and vary- 
ing displacements in the y-direction, relaxing the atoms 
(with the boundary layers contrained to be rigid and non- 
rotating) , and measuring the grain boundary energy per 
length. The range of displacements in the y-direction 
that we must search over is given by 

Di D2 mod Di^O 

A={ mm{\D2-[D2/Di\Di\, 

\D2 - \D2/Di^ Di\) D2 mod Di ^ 

(1) 

where we assume Di < D2. This is illustrated in figure[3l 
The grain boundary energy per length is defined as 



Egb — 



Etotal ~ Natorns * Ebulk 



(2) 



where Etotai is the total potential energy for the con- 
figuration of atoms (excluding the constrained atoms), 
Natoms is the number of unconstrained atoms, E^^ik is 
the energy of a single atom in the bulk, and L is the 
length of the grain boundary. An example of the results 
of such a search is shown in figurclH The regions with the 
same final displacement and energy correspond to basins 
of attraction around the finite number of final configu- 
rations of atoms. The minimized energy per unit length 
(equation found by this method is what we record as 
the grain boundary energy. 

We have also tried thermal annealing and have found 
that for certain geometries, the grain boundary migrates 
to form a jagged interface with segments of different grain 
boundaries that collectively have a lower total energy 
than the geometry given by the original set of tilt an- 
gles. The tendency of certain grain boundaries to corru- 
gate is also discussed by Ishida and Pumphrey^^iii. As 
a corrugated grain boundary is torn apart to measure 
the cohesive law, the corners form stress concentrations 
which weaken the grain boundary. Despite the fact that 
this configuration may be more natural, it is not what 
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FIG. 3; The Necessary Search Range for y- 

displacements. The minimum range we must search over 
to find the most natural configuration is equal to the minimum 
shift that gives an equivalent configuration of atoms along the 
surface. Assuming that the edges of the repeat cells line up at 
the bottom, the minimum shift to produce an equivalent con- 
figuration is done by lining up the top edge of the first repeat 
cell in the grain with a larger repeat (the grain on the right 
in the diagram) with the nearest edge of a repeat cell in the 
grain with a smaller repeat distance (the grain on the left). If 
D2 is equal to an integer number of repeats of Di, the mini- 
mum search range is equal to D\. YD2/ D\\ gives the number 
of repeats of the left grain that fit within a single repeat of 
the right grain and \D2/D{\ gives the number of repeats of 
the grain on the left that contain one repeat of the grain on 
the right. If D2/D1 is not an integer, the minimum search 
range is then min(lD2 - [D2/D1JD1I , \D2 - \D2/Di~\Di\) 



we intend to measure. We wish to measure the fracture 
toughness of all possible pairs of tilt angles, even if those 
pairs of tilt angles happen to be unstable. In a few cases, 
the minimization procedure described above produced a 
curved grain boundary. In these cases, we constrained 
the displacements in the y-direction such that a consis- 
tent pattern of flaws along a straight grain boundary is 
achieved. 

After finding the ideal y-displacement, we increment 
the strain by displacing the constrained layers of atoms 
in the a;-direction, away from the grain boundary. We 
relax the atoms and measure the force in the x-direction 
per unit length on the constrained layer of atoms. If 
the measurement of the stress drops abruptly during one 
strain step, the simulation restores the positions from a 
previous step, reduces the size of the strain increment, 
and proceeds. The result of such a simulation is shown 
in figure[T] The maximum stress in the stress strain curve 
is what we define as the peak stress. 



III. GRAIN BOUNDARY GEOMETRIES 

A. Lattice Symmetries and Tilt Angles 

For the triangular lattice, it is clear that we only need 
to explore tilt angles between 0° and 60°, but we can 
further reduce the space of grain boundary geometries 
to consider. Figure [S] represents the space of tilt an- 
gle pairs. Reflecting a point in this space through the 
62 = Oi line corresponds to swapping the two grain ori- 
entations which, as shown in figure [51 is equivalent to 
fiipping the grain boundary in both the vertical and hor- 
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FIG. 4: Finding the Energy Minimizing Grain Bound- 
ary Configuration. Each plateau in the top plot and each 
flat region in the bottom plot correspond to the basin of at- 
traction of a local minimum. The most natural grain bound- 
ary configuration corresponds to the global energy minimum. 



izontal directions, resulting in the same grain bound- 
ary. Reflecting through the (30,30) point takes (6*1, 6*2) 
to (60 — 6*1,60 — 62), reversing the sense of rotation of 
each grain. As shown in figure [SJ this is equivalent to 
flipping the grain boundary in the vertical direction. Re- 
flecting through the 62 = 60 — 9i line takes {01,62) to 
(60 — 6*2, 60 — 6*1), both switching the grains and reversing 
the senses of rotation. This is equivalent to flipping the 
grain boundary in the horizontal direction, also shown 
in figure [5l Thus, we only need to consider the pairs of 
tilt angles in the triangle enclosed by the 6*2 = line, the 
6*2 = 9i line, and the 6*2 = 60 — 6*1 line - the shaded region 
in figure [51 



B. Finding All Possible Geometries for Periodic 
Boundary Conditions 

In order to better simulate the grain boundary in the 
bulk as described in section |TT1 we would like to use peri- 
odic boundary conditions along the direction of the grain 
boundary (y-direction) . In order to simulate 2D grain 
boundaries in periodic boundary conditions, we need not 
only for both grain orientations to have finite repeat dis- 
tances, but also for the repeat distances to be commen- 
surate with one another. Let (a, b) be the lattice vector 
that is parallel to the edge of a 2D triangular lattice. 




FIG. 5: Finding the Minimum Necessary Range of 
Angles. For each point in the shaded region, there is an 
equivalent grain boundary in each of the other regions. 

Since the basis vectors arc at a 60° angle to one another, 
the repeat distance of a particular orientation of a 2D 
triangular lattice is given by 



D = y (a + -&)2 + (^6)2 ^ +p + ah. (3) 

The tilt angle is then given by 

a --if a + h/2 \ 

9 = sm . (4) 

VVa2 + &2 + a5y 

In order to expand the number of possible geometries, 
we have also considered geometries for which applying 
a small strain to each grain allows us to fit both grains 
inside the same periodic box. For any pair of orientations 
that have repeat distances that are not commensurate, 
we can find a continued fraction approximation to the 
ratios of their lengths, p/q « Di/D2- We can strain each 
grain equally into a box of size L = {l/2){pD2 + qDi), 
where the strain required is \L — pD2\/ L. For a strain of 
0.05%, the grain boundary would have to be 1,000 A long 
before the strain would alter the structure of the grain 
boundary. We find all pairs of tilt angles that correspond 
to commensurate or near commensurate grain boundaries 
by looping over pairs of surface vectors and comparing 
the repeat distances. 

Figure [S] shows all possible geometries that can be sim- 
ulated with a periodic length of 70 lattice constants or 
less, and a strain of 0.05% or less. The 9i — 6*2 line cor- 
responds to perfect crystals, while the 6*1 = 60 — 02 line 
corresponds to symmetric grain boundaries. There is a 
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FIG. 6: All Possible Geometries. The set of points 
above represent all 2D grain boundary geometries that can 
be simulated in a periodic box of 70 lattice constants or less, 
with a strain of 0.05% or less. There are gaps near per- 
fect crystals, symmetric grain boundaries, and high symmetry 
grain boundaries (discussed in section UlI C|l because creating 
a new, nearby geometry requires adding flaws at large dis- 
tances. Lines radiating from high symmetry geometries rep- 
resent flaws added at closer and closer distances as you move 
away from the point representing the original geometry. Lines 
of slope ±1 represent a constant misorientation between the 
two grains. 

gap near each of these hnes because creating very small 
angle grain boundaries, or a geometry very close to sym- 
metric grain boundaries, requires adding flaws that are 
separated by large distances. We see gaps near other 
high symmetry grain boundaries for the same reason. 
High symmetry grain boundaries are discussed in sec- 
tion IIII CI The lines radiating out from high symmetry 
grain boundaries represent adding single flaws to those 
high symmetry grain boundaries at larger and larger dis- 
tances, as you approach the high symmetry boundary. 
Each line represents a different type of flaw. The sym- 
metry about the di = 30 line is due to the fact that 
grains with tilt angles 6 and 60 — 6' have the same re- 
peat distance. Therefore, the grain boundaries given by 
the tilt angles (0i,02) and (60 — 61,62) have the same 
overall repeat distance, though these are different grain 
boundaries. 



C. High Symmetry Grain Boundaries 

Certain grain boundary geometries have particularly 
low repeat distances. These grain boundaries have spe- 
cial properties. They mark the center of cusps in the 
grain boundary energy and discontinuous increases in the 
fracture strength as a function of tilt angle. Tablc[T]show 
examples of high symmetry grain boundaries and figurc[7] 
shows examples of geometries in between, which consti- 
tute adding a single flaw per supercell repeat distance to 
the high symmetry grain boundary. Note how the added 
flaws constitute a compromise between the two high sym- 
metry geometries. Furthermore, a less high symmetry 
grain boundary can be repeated, and again have a sin- 



gle flaw added. Because of this hierarchical procedure 
for constructing lower-symmetry boundaries, the grain 
boundary energy and fracture strength as a function of 
tilt angle will have a self similar nature. 

TABLE I; High Symmetry Grain Boundary Geome- 
tries 
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The coincidence site lattice model (CSL) describes 
grain boundaries in terms of S, the inverse density of 
lattice sites that are shared by the two grain orienta- 
tions when rotated about a common lattice point^. In 
2D, CSL grain boundaries are necessarily commensurate 
and vice versa. For the triangular lattice, each common 
lattice point for a given pair of grain orientations has a 
repeat cell that is an equilateral triangle with one edge 
defined by the surface vector with length D. E is equal 
to the number of lattice points inside this cell and is 
given by the area of the cell divided by the area of one 
lattice triangle, S = D^. While the CSL formulation 
gives a simple method for finding commensurate grain 
boundaries, it is misleading to suggest that the actual 
coincidence of sites plays a physical role. In fact, shift- 
ing the grains in the j/-direction so as to minimize the 
grain boundary energy (described in section [Tl| generally 
causes the atoms to no longer coincide. Without this 
shift, the two free surfaces will have atomic planes that 
meet at the same point. The elastic energy is lowered by 
staggering the dislocations^^. Wolf has pointed out that 
for 3D, the CSL formulation of grain boundaries involves 
a redundant number of parametersi^. 

Bishop and Chalmers introduced the concept of "struc- 
tural units" - polygonal structures of atoms along high 
angle grain boundaries'^. We find that the patterns of 
"structural units" are more relevant to the physical prop- 
erties of the grain boundaries. When relaxed, the struc- 
tural units of grain boundaries in 2D, triangular lattices 
arc 5 atoms forming a regular pentagon or a pentagon 
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FIG. 7: Finding the Burger's Vector for Flaws Along 
Vicinal Grain Boundaries. Grain boundaries with tilt an- 
gles in the vicinity of high symmetry grain boundaries are 
comprised of the pattern of structural units of the high sym- 
metry grain boundary with an element from the neighboring 
high symmetry geometry. When defining the Burger's vector 
of a flaw added to a high symmetry grain boundary we are no 
longer comparing the flaw to the perfect lattice but to the re- 
peating pattern of structural units. Therefore, we must cross 
the grain boundary at equivalent portions of the structural 
unit in our reference loop and in our loop around the flaw. 
The loop on the bottom is the reference loop around the re- 
gion of the grain boundary without the added flaw. The loop 
on the top surrounds the added flaw, which we model as a 
partial grain boundary dislocation. The arrows indicate the 
segments needed to make the top loop match the bottom loop; 
for example, the arrow on the upper right indicates that the 
path around the defect is one atomic length shorter on that 
leg than the corresponding leg of the reference loop. The ar- 
rows sum to the Burger's vector of the added grain boundary 
dislocation. This grain boundary has a tilt angle of 44.3° 
and the added partial dislocation has a Burger's vector with 
a length of 0.495 lattice constants and a direction along the 
negative s-axis. 

that is slightly stretched. Each high symmetry grain 
boundary is comprised of a simple pattern of pentago- 
nal structural units (table |T| . When a high symmetry 
grain boundary is perturbed, an element of the pattern 
of structural units from the neighboring high symmetry 
grain boundary is introduced (figure [7]). As the tilt an- 
gles move closer to the neighboring high symmetry grain 
boundary, the flaws become closer together until they 
outnumber the original structure. The roles of flaw and 
original structure are then reversed. In this manner, com- 
binations of patterns of structural units can be used to 
build up any commensurate, high angle grain boundary. 

We find the structural units a useful way of concep- 
tualizing the different grain boundary geometries. We 
have not found the coincidence site lattice model useful in 
our investigations. In systematically constructing grain 
boundaries with small repeat distances (scction|TTTB]), wc 



have found the surface lattice vectors of the two sides of 
the grain boundary to be a useful description. We sug- 
gest that the structure of high angle grain boundaries can 
best be described in terms of a dislocation model. The 
"extra flaws" added to create vicinal grain boundaries 
can be described as partial dislocations. Figure [7] shows 
how to find the Burger's vectors for the added flaws by 
examining the flaw in the pattern of structural unitsj^ 
We will show in section IV Bl that this model gives excel- 
lent agreement in the stress field due to the added flaws, 
outside the background of the original fiaw structure. We 
will show in sections ITVl and IVl that the dislocation model 
provides the most powerful framework for understand- 
ing the geometry dependence of the properties of grain 
boundaries. 



IV. GRAIN BOUNDARY ENERGY 

We have measured the grain boundary energy and peak 
stress for all symmetric grain boundaries with repeat dis- 
tances under 20 lattice constants and asymmetric grain 
boundaries with repeat distances under 30 lattice con- 
stants. In order to explore the regions close to high sym- 
metry grain boundaries, we have added a few geometries 
with longer repeat distances, close to the high symmetry 
grain boundaries. 

The energy^ associated with the series of flaws that 
make up the grain boundary is deflncd in equation [21 
As stated earlier, and found by several earlier stud- 
jgg7|9ji6ji7ii8ii9i22 ^ cusps appear at high symmetry bound- 
aries. For small angle grain boundaries, the grain bound- 
ary energy has the form 

where fj, is the shear modulus, b is the Burger's vector, v 
is the Poisson ratio, and a is a factor that includes the 
core energy2£. One can now imagine the same scenario 
applied to high symmetry boundaries. We can take a 
high symmetry grain boundary and add or subtract flaws 
a distance d apart as shown in figure [5] By the same 
reasoning as used for the low angle grain boundaries, the 
energy near the high symmetry grain boundary will have 
the form 

where Eq is the energy of the high symmetry grain 
boundary which occurs at the angle a-nd a' incorpo- 
rates the core energy of the flaw within the pattern of 
flaws. 

The grain boundary energies for all of the symmetric 
grain boundaries that we have measured are shown in 
figure [51 Note that cusps occur at the angles listed in 
table m We are able to fit the data for symmetric grain 
boundaries to a function of the form 
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FIG. 8: Adding a Flaw to a High Symmetry Grain 

Boundary. The lighter dislocations represent flaws added a 
distance d apart, to an existing pattern of dislocations, shown 
in black, with a short repeat distance. The added flaws can 
also move, screen, or cancel the flaws that make up the high 
symmetry boundary. 
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COS 69i I log 1 TTa^ — TTa-r was chosen as a fitting function 

*l O I cos DP — cos D^i I ^ 

because it asymptotically gives a 6'log(l/6') shaped cusp 
in the near vicinity of 9i and because it has the correct 
symmetry: even mirror symmetry at 0° and 30° and an 
overall period of 60° . We use one term that is antisym- 
metric about 9i and one term that is symmetric about 9i 
so that we can fit the shape on either side of the cusp in- 
dependently. Wc do not expect the slope of the curve on 
either side of the cusp to be the same since the Burger's 
vectors of the additional flaws for the geometries on either 
side of the high symmetry geometry may differ (figure[7]). 

The 9i can be any angles that have the shortest re- 
peat distances, such as those given in table HI The 
curve in figure [H is the result of fitting equation [7] to 
the data shown in the same figure. We have used 
9i = (49.10°, 43.89°, 40.89°, 38.94°) since these angles 
have particularly prominent cusps and short Miller in- 
dices (given in table Three smooth sinusoidal terms 
were used in the final sum over cos 6 j9. 

The result is analogous to a devil's staircase, with a 
cusp singularity at each angle that corresponds to a spe- 
cial rational number. In principle, the energy has a loga- 
rithmic cusp at a dense set of points but (as is typical for 
devil's staircases) the high-order cusps rapidly diminish 



FIG. 9: Grain Boundary Energies for Symmetric 
Geometries. Cusps appear at high symmetry grain bound- 
aries (listed in table |l]) and have the same 6log0 shape as the 
energy of low angle grain boundaries. The line is the fit given 
by equation [T] Notice that hints of smaller cusps are visible 
in the data. 
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The first two terms fit the cusps at 0° and 30°, where 
the cusps are symmetrical about their respective center 
points. The next set of terms in the sum fit the cusps 
at high symmetry tilt angles. The function | cos 60 — 



TABLE II: Coefficients for Fitting Grain Boundary 
Energy to Eq. (T) 
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The results for the energy and peak stress of all of 
the grain boundaries we have simulated are shown in fig- 
ure [TOl As with the symmetric grain boundaries, we see 
energy cusps for grain boundaries with particularly low 
repeat distances. For example, one of the cusps in fig- 
ure [To] is centered around the grain boundary with tilt 
angles 24.79° and 3.00° which has a repeat distance of 
9.53 lattice constants. Figure [Til shows geometries in the 
neighborhood of this high symmetry boundary with a 
constant misorientation. We see a similar, 9\og9 shape 
to the cusp centered around the high symmetry geom- 
etry, suggesting that the same concept of adding flaws 
to the short-repeat distance grain boundary applies to 
asymmetric 2D grain boundaries along constant misori- 
entation lines. 
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FIG. 10: Grain Boundary Energies for Asymmet- 
ric Geometries. The surface and contour represents a 
smooth interpolation between the particular geometries that 
have been simulated, represented by the black dots in the 
6182 plane. Perfect crystals have zero grain boundary energy, 
shown by the groove along the 61 = 62 line. 
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FIG. 11: Asymmetric Grain Boundary Energies in 
Neighborhood of High Symmetry Geometry. The grain 
boundary geometries considered in the above plot each have 
a misorientation of 21.79° and are in the vicinity of the high 
symmetry grain boundary with tilt angles 24.79 and 3.00. The 
lines are fit to 9 log 9 for the points on either side of the cusp. 



V. FRACTURE STRENGTH 

A. Low Angle Symmetric Grain Boundaries 

The Frank conditions state that the total Burger's vec- 
tor for the dislocations making up a low angle grain 
boundary is equal to the difFcrence of the surface vec- 
tors that define the orientation of each grain. In order to 
guarantee that the grain boundary will have only one dis- 



location per repeat distance after it is relaxed, we must 
choose surface vectors that have a difference equal to a 
basis vector and have the same repeat distance. One pair 
of surface vectors is (2n -I- 1, —n) and (2n -I- 1, —n — 1). 
This gives a symmetric grain boundary with a single dis- 
location with tilt angles close to 0°, a Burger's vector 
equal to (0, 1) and a repeat distance of vorz^ + Sn + l. 
Our simulations show that such grain boundaries fail via 
intragranular fracture rather than intergranular fracture 
as shown in figure W]\ 




FIG. 12: Intragranular Fracture for Low Angle Grain 
Boundaries. The figure on the left shows a grain boundary 
with symmetric tilt angle 0.81° at strain. The same grain 
boundary is shown on the right with a strain of 3.125%. Sym- 
metric low angle grain boundaries centered around the 0° ori- 
entation (Miller indices (0, 1)) fail via intragranular fracture 
rather than intergranular fracture. 

The apex of the triangle 6*1 = 6*2 — 30° in figure [6] is a 
perfect crystal oriented at 30° , which is the end point of a 
series of symmetric tilt boundaries. At exactly 9i ~ 02 ^ 
30°, there is an abrupt jump up in fracture strength since 
the perfect crystal has no "extra dislocations" at which 
to nucleate fracture. Our simulations find a peak stress 
of 4.31 (Lennard-Jones Units) for the perfect crystal. 

The low angle grain boundaries near the 30° lattice 
orientation have surface vectors and (n, 1), repeat 

distances and total Burger's vector (—1, 1), 

which splits into two fiaws with Burger's vectors (0, 1) 
and (—1,0) shown in figure [T5] For wide enough sim- 
ulations, these dislocations glide in opposite directions 
until they are restricted by the constrained zones on ei- 
ther side. For narrower simulations, the dislocations do 
not glide but form nucleation points for grain boundary 
fracture causing an abrupt jump down in the peak stress 
compared to the peak stress of the perfect crystal. 

For narrow simulations of low angle grain boundaries 
in this region, we find that the peak stress has a parabolic 
dependence on angle, shown in figure UM We can explain 
this parabolic dependence as due to a partial screening of 
the external stress <7ext by the neighboring dislocations 
on the grain boundary. Assume that the dislocation has 
a critical stress for nucleating fracture equal to ac- The 
dislocation feels a stress due to its neighboring disloca- 
tions, each a distance d apart, in addition to the external, 
applied stress. The total stress felt by each dislocation 
can be written 

oo 

neigh, disloc. n—0 
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FIG. 13: Intergranular Fracture for Low Angle Grain 
Boundaries. The first figure sliows a grain boundary witli 
symmetric tilt angle of 30.96 at strain. The second figure 
shows the same grain boundary at 3.625% strain. In this 
wider simulation, the dislocations glide apart (thin diagonal 
white stripes). The third and fourth figures shows a narrower 
simulation of the same tilt angle where the intergranular frac- 
ture nucleates at each dislocation. 
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FIG. 14: Peak Stress vs. Tilt Angle For Low An- 
gle Grain Boundaries. The peak stress has a large jump 
downward as soon as a grain boundary deviates from 6 — 30° 
(perfect crystal). Here the peak stress for the perfect crystal, 
at the 30° orientation is 4.31 (Lennard- Jones units), which 
would require a vertical scale 11 times as big. The disloca- 
tions forming the boundary act as nucleation sites for frac- 
ture no matter how far apart they are. After this jump the 
peak stress has a parabolic dependence on angle for low angle 
grain boundaries that are constrained by width to fracture in 
qualitatively similar ways. Fracture nucleates exactly at the 
dislocation and the first non-vanishing term in the stress at 
this point due to neighboring dislocations goes as 1 /(P , where 
d is the distance between dislocations. 



where <7ext is the external stress. The n — 1 term is the 
Volterra solution given by 



(Tyy{x,y) 



lib y{3x^ + y^) 
27t{1-v) (x2-Fy2)2 

fib yjx"^ - y'^) 



(8) 
(9) 



27r(l - v) (a;2 -f-?/2) 



jn 

2\2 



(10) 



where the x-direction is the direction of the Burger's vec- 
tor. Since the Volterra solution is odd, the n = 1 term of 
the stress at each dislocation vanishes as we sum over the 
neighboring dislocations on either side. The first nonva- 
nishing term in equation [7] is the n = 2 term which has 
three contributions. (1) The n > 1 terms are the mul- 
tiple expansions of the stress field^^ as well as nonlinear 
terms. (2) The nonlinear term in strain field has the 
form du/dx * dujdx, giving a power law of l/r2. (3) 
Geometrical restrictions cause some grain boundaries to 
have flaws unequally spaced in the y-direction (though 
for the results given in figure 1141 we have only explored 
geometries with equally spaced flaws). The grain bound- 
aries geometries used in figure 1141 do have flaws that are 
not aligned perfectly in the x-direction. In each of these 
cases, shifting the dislocation constitutes adding a dislo- 
cation dipole (adding one positive and one negative, can- 
celing a dislocation and adding a new one) , and therefore 
is another contribution to the Ijr'^ term. The external 
stress needed to produce a stress equal to CTc at each flaw 
is then 



0"c 



02 

d2 



(11) 



where 02 combines the three contributions described 
above. This parabolic dependence on 6* — 30° is depicted 
on the left-hand side of figure fTSl 



B. High Angle Grain Boundaries 

Figure [15] shows the results of the peak stress mea- 
surements for high angle symmetric grain boundaries. 
At the same points for which we had cusps in energy, 
we have discontinuous increases in fracture strength. By 
drawing the same analogy between adding dislocations 
to perfect crystals and adding flaws to high symmetry 
boundaries as described in figure [8] we can understand 
the discontinuities in the fracture strength at high sym- 
metry grain boundaries and the angular dependence of 
fracture strength near the high symmetry geometries. 

For high angle grain boundaries, the added flaw is no 
longer the sole nucleation site for fracture and fracture 
does not necessarily nucleate in the core of the added 
flaw. The added dislocation creates a stress field given 
roughly by the Volterra solution (eauation fTU)) with a pos- 
itive stress on one side, negative stress on the other, and 
a singularity at the center shown in figure 1161 The stress 
field differs slightly from the Volterra solution because 
the elastic constants of the material at the grain bound- 
ary vary from those of the perfect crystal. The fracture 
nucleates along the boundary in the region where the 
stress due to the added flaw is positive. Because frac- 
ture does not nucleate at the center of the added flaw, 
the Volterra solution as summed over the neighboring, 
added flaws docs not cancel at the nucleation site. This 
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FIG. 15: Peak Stress vs. Tilt Angle For High Angle 
Grain Boundaries. The peak stress as a function of tilt 
angle is discontinuous everywhere, with higher values at spe- 
cial tilt angles representing high symmetry grain boundary 
geometries. The dependence of peak stress on angle near the 
high symmetry grain boundaries (parabolic near 30°, linear 
near 49.11°) depends on the structure of the additional flaws 
that make up the nearby geometries. 



leads to a linear law for fracture strength as a function 
of tilt angle, for grain boundaries near high symmetry 
geometries. 




FIG. 16: Stress Fields Due to Dislocations. The figure 
on the left shows the stress fields surrounding two dislocations 
according to formula [TO] The figure on the right shows the 
stress fields surrounding the added dislocations as calculated 
according to the virial definition of atomic stress^. 



Consider the grain boundaries with tilt angles ranging 
from 49.39° to 53.41°, which are close to the high symme- 
try grain botmdary at 49.11°. The additional flaws that 
characterize these grain boundaries have Burger's vectors 
equal to (— '\/3sin6', — sin6') and (— -v/Ssinfl, sin0), where 
Q is the symmetric tilt angle of the lattice. The norms 
are 2|sin0| and the angles are 210 and 150°. For these 
geometries, the xx component of the stress field (due to 
two dislocations a distance D w fe/2(6'o — 6)) along the 
y-axis is 

^ /i(2^/3^/(go - d) - 3sin(go - d)) sin(go -^^K 

2^(1 - v){y{B^ -6)- V3sin(0o - ^^)) 



FIG. 17: Peak Stress vs. Tilt Angle Near High Sym- 
metry Grain Boundaries. There is a discontinuity in peak 
stress at tilt angles close to high symmetry grain boundaries. 
The plot above shows the peak stress for the grain bound- 
ary with tilt angle 49.1 (described in table [I} and the nearby 
geometries. 



11(3 - 2V3y){0o - 6) 
27ry(l-i.)(y- V3) 



+ 0{{eo-9)f (13) 



where Oq is the tilt angle of the high symmetry grain 
boundary. We need to look at simulations in the fixed 
displacement (narrow width) regime in order to observe 
where fracture nucleates^i We find that for geometries 
with this pattern of flaws, fracture nucleates at the same 
distance above the added flaw. The external stress 
needed to nucleate fracture at a distance y from the 
added flaw, along the grain boundary is then approxi- 
mately 



^peak — 



(14) 



which explains the linear dependence on angle shown in 
figures and [T71 

The results for the peak stress of all grain bound- 
aries with repeat lengths less than 30 lattice constants 
are shown in figure [THl We also note spikes in peak 
stress for grain boundaries with particularly low repeat 
distances. Consider again the grain boundary with tilt 
angles 24.79° and 3.00° which has a repeat distance of 
9.53 lattice constants. Figure [TOl shows the peak stress 
of this grain boundary and grain boundaries in its neigh- 
borhood, along a line of constant misorientation. 



VI. CONCLUSION 

We have systematically explored the space of 2D grain 
boundaries, the patterns of structural units that make 
up short repeat distance grain boundaries, and the flaws 
in the patterns of structural units that make up vicinal 
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FIG. 18; Grain Boundary Peak Stresses for Asym- 
metric Geometries As with figure [TOl the surface and con- 
tours represent a smooth interpolation between particular ge- 
ometries. Due to the interpolation, each high stress peak is 
rendered as a finite-width cone. The peak stress of perfect 
crystals is much higher than that of grain boundaries. In or- 
der to use a smaller vertical scale, we have not included these 
data points which would otherwise show a ridge along the 
diagonal. 
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tilt angle dependence of peak stress near special grain 
boundaries. 

In principle, with enough computer time it would be 
possible to conduct a similar study for the 5D space of 
idealized 3D grain boundaries. We suspect that the per- 
fect crystal analogy could also explain the cusps in energy 
and discontinuous spikes in peak stress for special 3D 
grain boundaries. The geometry dependence of energy 
and peak stress surrounding these special grain bound- 
aries will be more complicated since the grain boundaries 
can have a pure tilt, pure twist, or a mixed type of geome- 
try. In real polycrystals, the fracture strength and energy 
are further complicated by impuritites at the interface, 
emission of dislocations during fracture, and more com- 
plex geometries such as triple junctions of grains, which 
are often the nucleation site for fracture. Systematic 
studies of such systems are infeasible and better suited 
by on the fly simulations of local regions of interest. 
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FIG. 19: Asymmectric Grain Boundary Peak Stresses 
In Neighborhood of High Symmetry Geometry. The 

grain boundary geometries plotted above are the same as 
those in figure 1111 along a line of constant misorientation in 
the vicinity of a high symmetry geometry. The lines show a 
possible linear fit to the points on either side of the low repeat 
distance grain boundary. 



grain boundaries. We have shown that the patterns of 
structural units are key to understanding the singularities 
in energy and peak stress at special grain boundaries by 
drawing an analogy to perfect crystals. We have used 
this insight to find a functional form for the energies of 
2D symmetric grain boundaries and to understand the 
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^® For cracks propagating through a grain boundary, the di- 
rection of the crack front in relation to the structure of the 
grain boundary may affect the fracture toughness, but this 
dependence is not included in cohesive laws. 
Another alternative is to use a layer of constrained atoms. 
Typically, in a finite element simulation of a subsection of 
a material, the nodes on the surfaces of the model are con- 
strained to not move in a direction perpendicular to the 
surface. This is known as "roUered" boundary conditions. 
We can imitate these boundary conditions in an atomistic 
simulation by constraining a layer of atoms at the surface 
or edge of the simulation to not move perpendicular to that 
surface or edge. This will partially ameiliorate the edge ef- 
fects, suppress the Poisson effect, and has the advantage 
of allowing us to simulate grain boundaries of any geome- 
try^. We have compared simulations of different sizes with 
roUered boundary conditions in the j/-direction to a simu- 
lation of the same geometry, but with periodic boundary 
conditions in the j/-direction. We see that if we use roUered 
boundary conditions, we need a simulation that is 48 times 
longer in the j/-direction in order to converge to within 2% 
of the results found by a simulation with periodic boundary 
conditions. Hence, for computational efficiency, we choose 
to use periodic boundary conditions in the j/-direction. 
The Frenkel-Kontorova model^^ describes a one dimen- 
sional chain of atoms, connected by springs, subject to an 
external, sinusoidal potential. If the relaxed length of the 
spring and the width of the potential wells are commen- 
surate a minimum energy state exists. For the incommen- 
surate case, there will be a pinned phase with many local 
minima and an unpinned phase, depending on the depth 
of the wells in relation to the stiffness of the springs. 
Since this Burger's vector may be a sum of lattice vectors 
from either side of the grain boundary, the total may not 
be a full translation vector of the triangular lattice. 
We use 'energy' for 'energy per unit length' throughout. 
For simulations of a few angstroms wide, the crack opens 
up along the grain boundary with a slow unzipping mech- 
anism, starting at particular flaws. For wider simulations, 
the crack snaps open, soon after a peak stress is reached. 
The slow unzipping of the narrow simulations is charac- 
teristic of fixed displacement boundary conditions - since 
the displacements of a region close to the interface is fixed, 
the opening of the crack is controlled, allowing us to see 
the details of how the crack opens. The snapping open, 
seen with a wider simulation, is characteristic of fixed force 
boundary conditions. In the wider simulations, the mate- 
rial on either side of the grain boundary acts as a spring, 
which effectively softens the boundary conditions, approx- 
imating fixed-force boundaries for long lengths. When the 
constrained region is far away, the interface snaps open 
once it reaches the maximum stress that it can sustain. 



